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ABSTRACT 

Diskoseismic c-modes in accretion discs have been invoked to explain low-frequency 
variabilities observed in black-hole X-ray binaries. These modes are trapped in the 
inner-most region of the disc and have frequencies much lower than the rotation fre- 
quency at the disc inner radius. We show that because the trapped waves can tunnel 
through the evanescent barrier to the corotational wave zone, the c-modes are damped 
due to wave absorption at the corotation resonance. We calculate the corotational 
damping rates of various c-modes using the WKB approximation. The damping rate 
varies widely depending on the mode frequency, the black hole spin parameter and 
the disc sound speed, and is generally much less than 10% of the mode frequency. 
A sufficiently strong excitation mechanism is needed to overcome this corotational 
damping and make the mode observable. 
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1 INTRODUCTION 

Diskoseismic oscillations of accretion discs around relativistic objects have been studied for over two decades (e.g., Kato & 
Fukue 1980; Okazaki et al. 1987; Nowak & Wagoner 1991, see Wagoner 1999, Kato 2001 for reviews), and have been used as 
models for the time variability and quasi-periodic oscillations (QPOs) in X-ray emissions from black-hole X-ray binaries. 

The c-modes (or so-called corrugation waves) were first proposed to explain low-frequency variabilities as their oscillation 
frequencies are lower than the associated g-modes and p-modes (see section 3). Kato (1983) and Okazaki & Kato (1985) showed 
the existence of one-armed (m = 1), low-frequency modes in nearly Keplerian (Newtonian) discs, while later work (Kato 1989; 
Silbergleit et al. 2001) demonstrated the presence of low-frequency c-modes in relativistic accretion discs, particularly the 
one-armed corrugation waves with a single node (n = 1; see section 3 below) in the vertical direction, which oscillate at 
(approximately) the Lense-Thirring precession frequency evaluated at the outer edge of the trapping region. 

In a recent paper (Lai & Tsang 2008), we studied the global corotational instability of non-axisymmetric p-modes (n = 0) 
in black hole accretion discs. The mode is trapped inside the corotation resonance radius rc (where the wave pattern rotation 
speed uj/m equals the disc rotation rate fi) and carries a negative energy. We showed that when the mode frequency ui is 
sufficiently high, positive wave energy is absorbed at the corotation resonance, leading to the growth of mode amplitude. The 
mode growth is further enhanced by wave transmission beyond the corotation barrier. Non-axisymmetric g-modes, on the 
other hand, may contain a corotation resonance in the wave zone. Kato (2003) and Li, Goodman & Narayan (2003) showed 
that such g-modes are heavily damped as the wave propagate through the corotation resonance (see also Zhang & Lai 2006). 

Diskoseismic c-modes are trapped in the inner most regions of black hole accretion discs. Although their primary wave 
zones are separated from the corotation resonance, the wave can tunnel through the evanescent barrier and propagate again 
around the corotation. In this paper we calculate the analytic damping rate of c-modes due to wave absorption at the corotation 
resonance. In section 2 we briefly review the basic properties of perturbations in a thin isothermal disc, and present the linear 
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perturbation equations. In section 3 we discuss the propagation regions associated with various diskoseismic modes, while in 
section 4 we demonstrate the effect of the corotation resonance on wave propagation. In section 5 the effect of the corotation 
on the c-mode is studied and the c-mode damping rates are calculated for different disc parameters. Section 6 contains our 
conclusion. 



2 BASIC SETUP AND EQUATIONS 

Consider a thin isothermal disc with the unperturbed velocity Uo = (0, rO., 0) in the cylindrical coordinates. The vertical 
density profile is given by (for small z <^r) 

Po(r,2) = ^^exp(-zV2^^'). (1) 
V ItxH 



Here E(r) is the (vertically integrated) surface density and H — Cs/Q.± is the vertical scale height, where Q.± is the vertical 
oscillation frequency and Ca is the isothermal sound. 

Perturbing the mass and momentum conservation equations gives 

^^5p + V -{poSu + vLoSp) = 0, (2) 

^(5u+(uo-V)<5u+((5u- V)uo = -VSh, (3) 

with the enthalpy perturbation 5h = 5P/ p — c^Sp/ p, where we assume that the perturbations are also isothermal. Assuming 
perturbations of the form 5P, 5u, 5p oc eicp{im(f) — itut), we have 



^6h + - — {rpoSur) -\ poSu^ + —(pSu^) = (4) 

—iCjSur — 2nSud, = —-^Sh (5) 
or 

—icjSu^ + ^y'^'"'- ~ ~" — (6) 
d 

—iCbSuz = ——5h (7) 
oz 

where u = ui — mil. Following Okazaki et al. (1987), we assume a ^-dependence of the perturbations such that Sh, 5ur,5u^ oc 
H„(z/_ff) where }ln(z/H) is the Hermite polynomial of order n. Then equation (|4]) reduces to 



tu!^dh-\ — T^irpoOUr) OoOus —oh — (8) 

d r or r luj 



Neglecting terms proportional to dH/dr ~ 0(l/r), and eliminating the velocity perturbations 5ur and Sutf, from equations 
©, dSl) and (O, we obtain (see eq. [29] in Zhang & Lai [2006]) 



—Sh - (—In—] —Sh + 
dr^ \dr rTiJ dr 



rCu dr f2E 



Sh = , (9) 



where D = — co^. This is our basic working equation. 

The approach above is Newtonian. More rigorous fully relativistic derivations of the dispersion relation have been given 
by Ipser (1994, 1996), Perez et al. (1997), and Silbergleit et al. (2001). Some aspects of the general relativistic effects can 
be incorporated into our analysis by using the Paczynski-Witta psuedo-Newtonian potential, which gives n < — For 
our purposes of estimating the c-mode damping rates, it suffices to employ equation Q but with the relevant fully general 
relativistic frequencies (e.g., Aliev & Gal'tsov 1981; Okazaki et al. 1987) 

- ^{^~^+'-^y\ (11) 

6 8a Sa^ ^ 



where the frequencies are in units of c? /GM, r in units of GM/(? , and a is the spin parameter of the black hole. 
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Figure 1. The propagation diagram for non-axisymmetric g-modes and c-modes. Note that the curves mQ — k and mQ + k meet at 
r = rjsco since K(rjgco) = due to general relativistic effects. The g-modes are trapped between the inner Lindblad resonance (where 
LU = mQ — k) and the outer Lindblad resonance (where lu = mQ + k), or around the peak of the mQ + k curve. 

3 PROPAGATION DIAGRAM AND C-MODES 

There are three possible critical resonant points in the disc: the Lindblad resonances (LRs) where D = 0, the vertical resonances 
(VR) where uj^ — nQ^, and the corotation resonance (CR) where uj = 0. Far from these critical points, the WKB dispersion 
relation [for 5h oc exp(i J kdr)] takes the form (Okazaki et al. 1987) 



The modes with n = have no vertical structure, and are referred to as p-modes - their stability properties are studied in 
Lai & Tsang (2008). We focus on the modes with n ^ 1 and m ^ 1 in this paper. Wave propagation is allowed in regions 
where < < nO]^ or uj'^ > nf25_ > (note that k < Q,j_ in GR). The former defines the g-mode propagation zone: 
the mode is trapped in the region where mQ — k. < uj < raO, + the latter leads to c-modes, for which the wave zone is 
specified by tJ < mil — ^/nQ.± (see Fig. 1) . Clearly, the c-modes exist only when mQ. — ^/nQ±_ > and wave refiection occurs 
at r = Tin = risco, the radius of the Inner-most Stable Circular Orbit. For Newtonian discs, since Q± = k, = Q the c-modes 
can only exist if m > ^/n; for relativistic discs, however, Q > Q±_ (for spinning black holes), we can obtain modes for m = y/n 
that may have very low frequencies. The "fundamental" c-mode, with m = n = 1 is of particular interest (Kato 1990), since 
it corresponds to the Lense-Thirring precession of the inner, tilted disc. Note that for extremely low mode frequencies, the 
terms of order 1/r previously ignored may become important, and care must be taken to obtain the real eigenfrequencies 
for trapped modes. Here we ignore these complications and refer to Silbergieit et al. (2001) for a more thorough relativistic 
analysis. 

4 WAVE ABSORPTION AT THE COROTATION RESONANCE 

The Lindblad resonances [D — 0) are apparent singularities of the master equation ((9)l [this can be seen easily by writing ((9)l 
as two coupled first-order differential equations], and no wave absorption occurs at the Lindblad resonances (e.g., Goldreich 
& Tremaine 1979; Li et al. 2003; Zhang & Lai 2006). The vertical resonances (where to'^ = nQ\^) act purely as turning points, 
and no wave absorption occurs there either. However, the corotation resonance must be treated more carefully (Kato 2003; 
Li et al. 2003; Zhang & Lai 2006). 
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Here we follow the analysis of Zhang & Lai (2006). Near the corotation (r = rc, where uj — mQ), equation Q can be 
written as 

since for a thin disc the sound speed Cs <^ ri} and the last term in © dominates the other terms. Defining x = {r — rc)/rc 
and expanding (|14p around a; = 0, we have 

where 

In equation (|15p . we have inserted a small imaginary part ie (with e > 0) in because we consider the response of the disc 
to a slowly growing perturbation. 

The two independent solutions to equation (|15|l are 

5h± = z^/^'z^"' = z'/^e*'" " (17) 



where v = ^ C — \ S> 1 and z = x + ie. The solution z^^'^z"' has a local wavenumber k = ci(!^ In 2)/dr = u/{rcx), with the 
group velocity Vg = duj/dk = —uj/k = —qrcX^io /v < (where we have assumed f2 oc r"'', with g > 0), thus it represents waves 
propagating toward small r. Similarly, the solution z^^'^z'^" has Hg > and represents waves propagating toward large r. 

As shown in Zhang & Lai (2006), waves with n ^ 1 can propagate into the corotation region and be absorbed there. 
Consider an incident wave propagating from the a; < (or r < Tc) region toward a; = 0, with the amplitude (up to a constant 
factor) 

5h{x <0) — z^^'^e'^"'^^'^ — ie"" {—xy'^'^e~^"'^'^^~^\ (incident wave). (18) 
The transmitted wave is simply 

Sh[x > 0) = a;^''^e~"''"^, (transmitted wave), (19) 

and there is no reflection. Thus the amplitude of the transmitted wave is decreased by a factor of e"'^" (Zhang & Lai 2006). 
Similarly, a wave incident from the r > rc region toward rc encounters the same attenuation. Since v ^ 1 for thin discs. We 
readily conclude that all waves incident upon the corotation will be absorbed (Kato 2003; Li et al. 2003; see Zhang & Lai 
2006 and Lai & Zhang 2008 for applications of this result to the problem of wave excitation by external forces). 



5 COROTATIONAL DAMPING OF C-MODES 

The result of section 4 shows that waves propagating through the corotation are heavily damped. Since g-modes trapped 
between the inner and outer Lindblad resonances must cross the corotation, they are damped very quickly by corotation 
absorption, as shown by Kato (2003) and Li et al. (2003). Only higher frequency g-modes which are trapped around the peak 
of mO. + K (see Fig. 1) can avoid such corotational damping; these modes have frequencies ui ~ mr2(risco). 

In this section we calculate the damping rate of the c-mode. The damping mechanism is illustrated in Fig. 2, where we 
also plot the effective potential: 

Kff 2-2 + — + —^\l-^^-Fr)- 20 

ciuj^ r^ ruj \dr D / 

Based on the result of section 4, we will adopt the approximation that waves transmitted through the barrier between the 
inner vertical resonance (IVR) and the inner Lindblad resonance (ILR) are absorbed at the corotation. We will also assume 
that the inner disc boundary is completely reflective (see section 6). 

5.1 Reflection Coefficient 

We first calculate the reflection coefficient TZ when a wave in region I (see Fig. 2) propagates outward and is reflected back at 
the IVR. 

From the dispersion relation (|13p . the group velocity of the wave is given by 

dw kc? 
Ik " 



9 = -n-- -r T- (21) 



Q 1 - {K/Cj)^(nD.l/Cj^) 
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Figure 2. The wave propagation diagram (upper panel) and effective potential (lower panel) for c-modes. The modes are trapped in 
region I (between rjn = rjgco ^'^id the inner vertical resonance rjva), but can tunnel through an evanescent barrier (region II, between 
the IVR and the inner Lindblad resonance ^ilr) and propagate into the corotation (rc), where they are absorbed. 



The relative sign of Vg and the phase velocity Vp — uj/k is important. In region III (riLR < r < Vc), Vg and Vp have the same 
sign, thus waves propagating outwards correspond to fc > 0. In region I (r < rivft), Vg and Vp have opposite signs, so that 
waves propagating outwards have fc < and waves propagating inwards have fc > 0. 

As shown in the previous section we can, to good approximation, assume that waves propagating into the corotation are 
completely damped. Thus, only an outward-propagating wave exists in region III (see Fig. 2), with the wave amplitude (up 
to a constant prefactor) given by 

5h = Aexp(if kdr+jj, (22) 



where fc > is given by equation (|13p and and A = ^ D jrkll is the WKB amplitude. The connection formulae for the ILR 
(Tsang & Lai 2008) are 



\A exp (- \l^^'^\k\ dr) for r < hlr 

A exp f jfc| dr) for r < hlr 

^JJ kdr + 7r/4^ for r ^ hlr 



[ A exp I 



These then give for the evanescent zone (region II in Fig. 2); 



ILR \ / /"^'ILR. 



5h ~ ^ exp ( — / |fc| dr I + exp ( / |fc|dr 
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Figure 3. The reflection coefficient, 'R., of waves incident upon the inner vertical resonance (IVR), as a function of the (real) wave 
frequency LOr (in units of mf2in, where $7;^ is the disc rotation rate at r = rjn), for various black hole spin parameters (a) and sound 
speeds = PrVl\^. The upper panels are for m = 2, n = 1 and the lower panels for m = n = 1. 



— exp(-eii) exp 



k\dr \ + iAexp(+On) exp 



\k\dr 



where 

/■nLR 
eii= / \k\dr. 



The connection formulae at the IVR can be similarly derived; they are0 

5h^ ^ / 5^ (- /nval'^l '^') » '"'^^ 

\ A sin (J^'^'^ kdr + n/A) for r < nvR 



[ Acos(/;"'^fcdr + 7r/4) 



for r > rivR 
for r < rivR 



(25) 
(26) 

(27) 
(28) 



^ Note that around the IVR, the differential equation is matched by the Airy functions, rather than by the Airy function derivatives as 
in the case of the ILR. 
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(Abramowitz & Stegun 1964). Thus, we find that for r < nvR (region I in Fig. 2), 



oh ~ — e cos 



kdr + ■ 



+ i2Ae^''^^ sin 



kdr + 



Expressing this in terms of traveling waves and defining y — kdr — 7r/4, we have 



5h ~ iA 



-iy ^g+®ii 



1 -8: 
4" 



(29) 



(30) 



where the first term (oc e '^) corresponds to the incident (outgoing) wave and the second term (oc e*^) the inward going wave 
reflected from the IVR. Thus the reflection coefficient is 



7^ : 



e©ii + ie-ei: 



and the transmission coefficient (through region II) is 
T 



(31) 



(32) 



We can clearly see that \TZ\'^ < 1. 



5.2 Trapped C-modes and their Damping Rates 



Assuming a reflective boundary exists at r = r^^ < '"ivr, we can develop trapped c-modes in the inner disc, between r^n and 
the IVR. To illustrate this we consider a simple boundary condition at r = r^^'- 



Sh{nn) = 0. 

From section 5.1, the wave in region I can be written as 



5h — A exp{—iy) +TZ A exTp{iy), with y 



k dr — ir / 4. 



(33) 



(34) 



where k = kr + iki is complex. Applying the boundary condition (|33p to equation (|34p yields the eigenvalue condition: 

/■nvR 

exp(2ie) = -^j7^|, with Q= I kdr = Qr + iQi, (35) 



where Or and Qi are real. The real eigen frequency is given by 

''IVR 



e,. = 



kr dr — 



'■IVH ^(^2-ii2)(nf7i-ii2) , 3^ 
— dr = nn + — , 

CslUJr] 4 



(36) 



where = 0, 1, 2, ■ ■ • is an integer and ujr = uJr — mQ < in the trapping region between rin and rivR. The imaginary part of 
the frequency uji is determined by \TZ\ = cxp(— 20i), or 



tanhOi = — 



I7^| 



|7^| + 1 

Note that Qi = f'^" fci dr < 1 and 



dk 

' dio 



kr Cs 



UK Q, 



We then obtain 



u}i ~ -- exp(-26ii) 



1 - 



dr 

Cs 



(37) 



(38) 



(39) 



Note that in the trapping region (between rin and tivr), nAt^fi^/u;'' < 1. Thus the mode is always damped (uji < 0). 



5.3 Numerical Results 

Figure 3 depicts the reflection coefficients for waves impinging upon the inner vertical resonance (IVR) from r < rivR as a 
function of the wave frequency. We consider both the m = 2, n = 1 and the m = 1, n = 1 modes. The real frequency ranges 
from to (mQ. — \Ai^i)ki„- Waves with higher frequencies are "protected" by a larger potential barrier in the evanescent 
zone (larger Qu) and have \TZ\^ closer to uinty. 
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Figure 4. The real and imaginary frequencies for the primary trapped c-mode with the (5h(rin) = boundary condition and rin = risco 
for various sound speeds versus the black hole spin parameter. Note that there is no mode for the one-armed corrugation wave (the right 
panel) for sound speed Cs = O.lrQin. 



Figures 4-5 show the real and imaginary c-mode frequencies, computed using the WKB expressions derived in section 
5.2. Waves with higher frequencies have smaller damping rate |tJi|, consistent with the reflectivity results shown in Fig. 3. 
The damping rate is also smaller for cooler (smaller (5) discs. For the m > n case (e.g. Fig. 4, left panel), the dimensionless 
damping rate |tJi|/m57in increases with increasing a (while a;r/mSlin remans approximately constant), because the width of the 
region between tivr and riLR is smaller for larger a. For the m = n = 1 case (Fig. 4, right panel) the range of possible c-mode 
frequencies is bounded from above by the Lense-Thirring precession frequency at the inner boundary. Thus tj^ — *■ when 
a ^ 0. This leads to larger \LOi/uJr\ for small a since the potential barrier in the evanescent zone decreases with decreasing 
LJr/mQ.iYi. Note that although we use the relativistic frequencies (|10|) - (|12p our calculations not fully relativistic, and the real 
frequencies calculated shown in Fig. 4, particularly for the m = n = 1 mode with large black hole spin parameter a, are 
correct only in orders of magnitude [Fully relativistic calculation of the real frequencies of these trapped modes was done by 
Silbergleit et al. (2001).] Also note that adopting different inner disc boundary conditions would lead to different real mode 
eigenfrequencies than presented in Fig. 4. 

Finally, we note that equation (I39p is valid only when there is no loss of wave energy at the inner disc boundary rin. In the 
presence of the rapid radial inflow at the ISCO, we would expect additional mode damping due to the leakage of waves into 
the plunging region of the disc (see Lai & Tsang 2008 for a study of such leakage for p- modes). Alternatively, a sufficiently 
strong excitation mechanism is needed to overcome the corotational damping and make the c-modes grow. 



6 CONCLUSION 

In this paper, we have shown that diskoseismic c-modes suffer corotational damping due to wave absorption at the corotation 
resonance. These modes are trapped between the inner disc edge and the inner vertical resonance (where uj — mQ = —^/nfl±, 
with m,n ^ 1), but can tunnel through the evanescent zone and leak out to the corotation zone where wave absorption occurs. 
The mode damping rates are generally much smaller than the mode frequencies, and depend sensitively on the disc sound 
speed and the black hole spin parameters. 

With this paper, we now have in hand a complete picture of how the corotation resonance affects various diskoseismic 
modes in black-hole accretion discs, at least in the linear regime: Non-axisymmetric g-modes are heavily damped at the 
corotation resonance (Kato 2003; Li et al. 2003), while p- modes (inertial-acoustic modes) can be overstable due to the 
corotational wave absorption (see Lai & Tsang 2008 and references therein). The corotational damping rates of c-modes are 
much smaller than those of g-modes. 

Diskoseismic c-modes have been invoked to explain low-frequency variabilities in black-hole X-ray binaries (van der Klis 
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Figure 5. The damping rate \LOi\ for various sound speeds and black hole spin parameters for c-modes with purely reflective boundary 
conditions at ri„ as a function of the real mode frequencies u!r. The upper panels show the results for the m = 2, n = 1 modes and the 
lower panels show the m = n = 1 modes. The sound speeds are taken to be constant, Cg = f3rQi^. 



2006; Remillard & McClintock 2006). Fu & Lai (2008) showed that the basic properties of c-modes are largely unaffected by 
the disc magnetic fields and thus these modes are present in real discs. The results presented in this paper show that in order 
for the c-modes to be observable, a sufficiently strong excitation mechanism is needed to overcome the corotational damping. 
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